Three-Component Fermi Gas in a one-dimensional Optical Lattice 
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We investigate the effect of the anisotropy between the s-wave scattering lengths of a three- 
component atomic Fermi gas loaded into a one-dimensional optical lattice. We find four different 
phases which support trionic instabilities made of bound states of three fermions. These phases 
distinguish themselves by the relative phases between the 2/cf atomic density waves fluctuations of 
the three species. At small enough densities and strong anisotropies we give further evidences for a 
decoupling and the stabilization of more conventional BCS phases. Finally our results are discussed 
in light of a recent experiment on 6 Li atoms. 

PACS numbers: 03.75.Mm, 71.10.Pm, 71.10.Fd, 



Ultracold multicomponent atomic Fermi gases have re- 
cently attracted much interest [1[ . In particular the exis- 
tence of several internal degrees of freedom might stabi- 
lize some exotic phases. In this respect recent theoretical 
investigations strongly support the formation of a molec- 
ular state made of bound states of N atoms. For instance 
quartet (N = 4) and trionic (N = 3) states have been 
predicted in both three and one dimensions in the context 
of cold atoms systems @, S 1, |, 1, 0, i, ©, How- 
ever, these first studies assumed at least an SU(2) sym- 
metry and even an SU(iV) symmetry between the species 
which may not describe accurately the experimental sit- 
uation at non-zero magnetic field. Indeed in a recent 
experiment, where a stable N — 3 component mixture 
of atoms in three different hyperfine states of 6 Li has 
been stabilized at small magnetic field llj, the s-wave 
scattering lengths a mn between the three species exhibit 
strong anisotropic behavior as a function of the exter- 
nal magnetic field. In view of the promising perspective 
to observe trionic bound states in a near future, a care- 
ful study of the generic asymmetry between the species 
is clearly most wanted. It is the purpose of this work 
to do so. To this end we will study a three-component 
fermionic gas with equal densities, pi. 2, 3 = p, loaded into 
a one-dimensional (ID) optical lattice of wavelength A 
and transverse size a±. Away from resonance and when 
the 3D scattering lengths \a mn \ <C (A, a±), the system is 
described with a Hubbard-like model with contact inter- 
actions 1121 : 
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where cj n is the creation operator for a fermionic atom of 

color n — (1, 2, 3) at site i and pi >n — c\ n Ci tn is the local 
density of the atomic specy n. The Hamiltonian ([I]) is 
an anisotropic deformation of the U(3) Hubbard model, 
obtained when U mn — U, whose phase diagram has been 



recently elucidated [6j. In this case, for an attractive 
interaction U < 0, a spectral gap opens for the SU(3) 
spin degrees of freedom and one- and two-particle exci- 
tations are gapped for incommensurate density p. The 
dominant fluctuations consist into gapless Atomic Den- 
sity Waves (ADW) and SU(3)-singlet trionic excitations 
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. When U 12 £ U 23 + U»i, the 
continuous symmetry of ([I]) is strongly reduced to U(l) 3 
and the resulting anisotropy has dramatic consequences. 
Indeed, on top of the previous symmetrical phase, we 
find by means of combined low-energy and density ma- 
trix renormalization group (DMRG) approaches [H El 
that there exists for incommensurate density p three dif- 
ferent ADW phases supporting trionic instabilities and 
even decoupled BCS phases. 

The (U, V) model. Let us first start with the simplest 
symmetry breaking pattern, U(3) — * U(2) x U(l), when 
two species, say 1 and 2, play an equivalent role. In this 
case U12 = U, U23 = E/31 = V and |T]) may be viewed as a 
two-component fermionic Hubbard model with coupling 
U between the species (1, 2) which interacts with a third 
specy 3 with coupling V. As it will be discussed later, this 
model captures the essential features of the generic case. 
In the weak-coupling limit, its low-energy effective theory 
can be expressed in terms of the collective fluctuations 
of the densities of the three species by the bosonization 
approach 13|. Introducing three bosonic fields <j> n {x), 
the density operators for each specy read as follows: 

P + W^) _ J_an[2k F x + V^4> n (x)], (2) 
a y 7T Tra 

where x = ia, a = A/2 is the optical lattice spacing, and 
hp = 2np/X is the Fermi wave-vector. The second and 
last terms of Eq. J2J) describe respectively the uniform 
and 2kp fluctuations of the density operator of specy n = 
1,2,3. In our problem the interaction is best expressed 
in terms of the collective fluctuations of the total den- 
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sity, described by a bosonic field <&o = En=i^i>)/^ 
and of the relative density, described by a two-component 
bosonic field $ = ($y , $jj where $y = (0i -<f> 2 )j\f2 and 
$x = (^l + 4>2 — 203)/Vw. In terms of these variables 
the effective low-energy Hamiltonian of the (U, V) model 
splits into three parts, 7i = Ho +H S + 7i m ix, where: 



A 



(3) 



is the Hamiltonian of a Luttinger Liquid (LL) describing 
the low-energy properties of the total density fluctua- 
tions. In Eq. ([3]), 60 is the dual field to $0, «o = vp/K 
denotes the density velocity (vf = 2tasin(fci?a) being the 
Fermi velocity), and K = (1 + 2{U + 2V)o/3t:vf)~ 1/2 is 
the Luttinger parameter. The Hamiltonian TL S accounts 
for the remaining (spin) degrees of freedom and reads: 

f*=ll.-L 

2g± 



cos V27T$|| cos V67r$^ — ^ cos v 871-$ 1 1, (4) 

with Ay = g\\ = -Ua/2%, \j_ = (U - W)a/6n and 
g± = —Va/2-K. Finally 7i m ; x couples spin and den- 
sity fluctuations with 7i m i X = A m i X d x $od x $± where 
A mix = y/2(U - V)a/3ir. When U = V, i.e. A mix = 0, 
the spin and density fluctuations separate at low en- 
ergy, and model (0| is the bosonized version of the SU(3) 
Gross-Neveu (GN) model studied in Ref. @. In all other 
cases, A m ; x 7^ 0, and the spin and total density de- 
grees of freedom do not decouple, due to the anisotropy, 
even though we are considering incommensurate densi- 
ties. However, as we will see, at weak-enough couplings, 
i.e. when | A m i x /27rw^| -C 1, thanks to the opening of a 
spectral gap for the spin degrees of freedom, the spin- 
density coupling Hmix has little effect and can be safely 
neglected. In this regime the low-energy properties of 
the (U,V) model are captured by those of H s that can 
be elucidated by means of a one-loop Renormalization 
Group (RG) approach. For generic values of the cou- 
plings (U,V) we find that (A M ,<7 M ), fx — (||,-L), flow to 
strong couplings and the three species are strongly cor- 
related. In the strong-coupling regime, the bosonic fields 
$(x) get locked and a spin-gap opens. We further dis- 
tinguish between two phases, Aq and A^, depending on 
the sign of V. The Aq phase is obtained for V < 
and = (0,0) whereas the A n phase is stabilized 

for V > with = (t/tt/2,0). In both phases 

the low-energy spectrum is an adiabatic deformation of 
that of the SU(3) GN model and consists into three 
kinks (and anti-kinks) \u) n ), n — (1,2,3) [l5| ■ Under 
the SU(2) group acting on the species (1,2), these three 
kinks decompose into a doublet (|^i), l^a)) and a singlet 
\i1J3) with masses and velocities (my, «m) and (m±,v±) 
respectively. Though their wave functions are different 



in the two phases, they are labelled by the same quan- 
tum numbers as those of the original lattice fermions 
c] . We thus find that the one- and two-particle excita- 
tions are fully gapped in Aq^ phases. As a consequence 
the equal-time Green functions, G n (x) = (ct n Cj+ Xin ), 
are short ranged with Gu 2 \(x) ~ sm(kFx)e~ m ^ v ^ x \ 
Gs(x) ~ sin (kFx)e~ mxv - L ' x >. Furthermore defining 

Pnm{x) = (Pl nm Pi+x,nm) with P\ nm = cj „cj m , WC find: 

P l2 {x) ~ Q~m^jL\x\ and p 31(2) ( x ) ~ e-™\\n\*\, so that 
neither the Aq nor the A w phase support BCS pairing 
instabilities. The dominant fluctuations rather consist 
into 2fcp ADW with correlations N nm (x) = (pi t nPi+x.m) 
and trionic excitations made of three fermions. 

Atomic density waves and trions. In Aq,tt phases, upon 
integrating out the spin degrees of freedom, local density 
operators simplify as: 

p d x <$> (x) 



/3tt 



■ A„ sin [2k F x + V47r/3$o(z)], (5) 



where the amplitudes Ai = A 2 = Ay and A3 = Aj_ are 
non-universal functions of the couplings (U, V) and are in 
general different. We thus find in both phases a power- 
law decay for the ADW equal-time correlations functions: 
N nm (x) ~ p 2 + A„A TO cos(2kFx)\x\~ 2K / 3 . However the 
two phases Ao and A n distinguish themselves by the rel- 
ative sign of the amplitudes A„. Indeed, we find that in 
the Ao phase AyAj_ > and consequently that the 2k^? 
ADW of the species (1,2) are in phase with that of the 
specy 3. In contrast, in the A n phase, we have Ay Aj_ < 
and the 2kp ADW of the species (1,2) are out of phase 
from that of the specy 3. On top of these ADWs, Ao^ 
phases support trionic excitations made of three fermions 
with binding energy Eb ~ m±Vj_. These excitations can 
also be distinguished in Ao.n phases but in a weaker 
sense. In Ao the dominant trions are characterized by the 
equal-time correlation function To(x) = (Tj ^To^+x) ~ 
To sin(fci?x)|a;| _ ^ +9 / K -'/ 6 which is quasi-long ranged. In 
A n the trionic wave function with maximal k F amplitude 
is obtained when two atoms (1, 2) at one lattice site i bind 
antisymmetricaly with the third specy 3 at two neighbor- 
ings sites i - 1 and % + 1: = cl A^l-i ,3 ~ c l+i, 3 )- 
Its equal-time correlation function is given by T w (x) = 
(Tl ti T„ ti+x ) ~ T^sm(k F x)\x\- ( - K+9 ^ K ^ e so that both 
symmetric and antisymmetric trionic correlation func- 
tions always display a power-law decay and only their 
amplitudes depend on phases: \Tq\ > \T W \ in Ao and 
\Ttt\ > I To I in A n . The key quantity that distinguishes 
between Ao and A^ phases is thus the relative sign of 
the 2kp amplitudes Ay, Aj_ of the local ADWs ©. In 
this respect, when going from the A* to the .Ao phase, a 
quantum phase transition (QPT) takes place on the criti- 
cal line V = where Ay and Aj_ vanish and change their 
relative sign. There are two different QPT depending on 
the sign of U. In the type-I transition with U > 0, all 
degrees of freedom become massless at the transition and 



3 




X X 



FIG. 1: (Color online) DMRG results for (U/t,V/t) = 
(—4, —2) and p — 5/12 in the Ao phase. Both one-particle 
Green functions G„ and BCS pairing correlations Pi 2 are 
short range, while trionic correlations decay algebraically. 
Note that symmetric trions dominate with |Tb| > |Pn-| and 
local densities of all species rii(x) are in-phase. 



the critical theory consists of three decoupled LLs. In the 
type-II transition for U < 0, a QPT occurs in the two- 
component LL universality class where my =/= and only 
m± vanishes. In this case, the specy 3 decouples from 
the two others which form well defined BCS pairs with 
quasi-long range pairing correlations P\2(x) ~ a. 
being some non-universal exponent. 

Strong-couplings and Trionic-BCS transition. So far 
we have neglected the spin-density coupling H m ix- At 
weak couplings, when |A m i x |/27njF <C 1, we find that the 
only effect of 7i m i x consists into a small renormalization 
of the low-energy parameters and do not modify qualita- 
tively the two-phase structure discussed above. At larger 
couplings, when |A m ; x |/27r7j^ 3> 1, the structure of the 
7i mix term strongly suggests that it may be responsible 
for a decoupling between the pair (1,2) and the specy 3 
leading, on top of Ao t7r phases, to two additional phases: 
a BCS phase where atoms (1,2) bind into pairs and even 
a fully gapless phase of three decoupled LLs. In the limit 
of large attractive \U\/t 3> 1 and repulsive V/t > 0, a 
trionic-BCS QPT occurs from an An phase to a decou- 
pled BCS phase in the (1,2) channel at small enough 
densities [16l ] . Apart from this case, the question of how 
do the four phases, Aq, An, BCS and LLs, interpolate 
in the strong coupling or low density regime is a difficult 
problem which requires a thorough numerical approach 
like DMRG calculations. 

Numerical simulations. In order to check the above 
theoretical predictions, we have performed extensive 
DMRG calculations for various densities 1/12 < p < 
5/12 and couplings —4 < U/t, V/t < 4. Simulations are 
done on open chains (up to 144 sites) keeping up to 1600 
states. The complete phase diagram will be published 
elsewhere [lij and we only report here our main find- 
ings. At sufficiently large densities and weak anisotropies 
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FIG. 2: (Color online) Same as Fig. H for (U/t, V/t) = (-4, 2) 
and p = 5/12 in the An phase. The only difference in that 
case is that antisymmetric trions dominate with \T n \ > \To\ 
and local densities n\ and ni are out of phase with 713. 



the DMRG results strongly support the two phase struc- 
ture, A® and An, predicted by the weak-coupling ap- 
proach. As an example Fig. 1 and 2 show our results 
for G n (x), P nm (x), Tq >1z (x), as well as the local density 
profiles n n (x) = (p% n ) f° r a density p = 5/12 and typical 
values of the couplings in the Ao and An phases. At small 
densities and larger anisotropies we observe a strong ten- 
dency toward decoupling. For example, by lowering the 
density at fixed couplings (U/t, V/t) = (—4,4), we find a 
QPT toward a decoupled BCS phase in the (1, 2) channel 
at densities p < p c ~ 1/4 [l7j |. 

General asymmetric model. We are now in a posi- 
tion to discuss the general case where Uu ^ U23 ^ f/31. 
The resulting phase diagram in the parameters space is 
rich and complex and will be presented in details else- 
where [l6j|. It can be shown that at large length scales, 
the low-energy theory is then equivalent to that of an 
effective (U, V) model. Since there are three inequiva- 
lent ways to define such a model, we find that, on top of 
the Aq phase, three inequivalent An(n,m) phases can be 
stabilized. The properties of each of these phases follow 
from those discussed above for the case (77,777) = (1,2) 
by a suitable permutation of the indices in the correla- 
tion functions. At large couplings and/or small densities, 
the system decouples and three BCS(t7, to) phases can be 
stabilized as well as a fully gapless decoupled LL phase. 

Experimental realization. A stable mixture made of 
a balanced population of three hypcrfine states of 6 Li 
atoms, \F,m F ) = |1) = |l/2, 1/2), |2) = |l/2, —1/2), and 
1 3) = 1 3/2, —3/2), has been stabilized recently in an op- 
tical dipole trap [III [Lsj . One may in principle further 
load the atoms in a 3D optical lattice with potential: 
V(x,y,z) — s±Efilsin 2 (kx) + sin 2 (fcy)] + s^Ensm 2 (kz) 
where sj_,|| = V q ±^\/Er, Er = h 2 k 2 /2M being the recoil 
energy. A ID optical lattice in the z direction would then 
be further stabilized by increasing the lattice potential to 
a high enough value s± 3> st\ and s± 3> 1. Neglecting 
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(n, m) is macroscopic and one finds that in the limit of 



large sample size L, N n>m ^ whereas N„ 



0. 
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FIG. 3: (Color online) Effective Hubbard parameters U nm as 
a function of magnetic field. The cross indicates the critical 
field B c between Ao and At, (or BCS) phases. 



the harmonic potential and for small enough scattering 
lengths a mn , the low-energy physics of such a system is 
captured by the fully anisotropic Hubbard model ([1]) [l2j 



with parameters U nrn = ^8/%E R (sj_S||) 1/4 aid, mn /a± 
and t = A/y/7rE R sjj^e -2 ^ where aid, mn = a mn /(l - 
(C/V2)(a mn /a±)) is the effective ID scattering length, 
a± = X/2ns , 1/ ' 4 the transverse confinement length and 
C = 1.4603 [ill- We show in Fig. [3] the dependence of 
the ratio U mn /t as a function of the external magnetic 
field B for typical optical lattice parameters A = lpva, 
s± = 20 and s\\ = 4. Using the one- loop RG approach 
discussed above and large scale DMRG calculations, we 
find the following phase diagram which is depicted in 
Fig. [3] An Ao phase with symmetric trions is stabilized 
independently of the density for magnetic fields B < B c . 
Above B c and at large enough densities p an ^4.^(2,3) 
phase emerges. The latter phase is unstable toward de- 
coupling when decreasing the density below p < 1/3. In 
the decoupled phase a BCS instability occurs with pairs 
of atoms in states 2 and 3, the specy 1 being decoupled. 
The critical field is estimated with the help of RG equa- 
tions to be B c ~ 563G, a value which is consistent with 
our numerical data. The numerical values of the trionic 
binding energy strongly depend on the phases. In Aq 
they are mostly independent of the density and only de- 
pend on B. For example, we find trionic binding energies 
E b /k B ~ 2600nK for B = 320G and E b /k B ~ lOOnK 
for B = 553G at all densities. In the A Tr (2,3) phase 
(i.e. B > B c and p > 1/3), we find that the trionic 
binding energies are small (typically Eb/kB < 30nK). 
In the decoupled case (i.e. p = 1/6 and B > B c ), we 
estimate the BCS gap to be of the order lOOnK. The 
different phases discussed above may be probed in ex- 
periments l^i 13] by measuring, with absorption imag- 
ing and via a series of magnetic field ramps, the average 
numbers of paired atoms (nm) relative to the non inter- 
acting theory: N n . m = 1/L dx[(p n (x)p m (x)} -p 2 ]. In 
a decoupled BCS phase with pairs in the (n, to) chan- 
nel and decoupled specy p, the number of bound pairs 



In both trionic phases all atoms are bound into pairs and 
N m ,n ^ 0, N m . p ^ and N p . n ^ 0. Though in the Aq 
phase all iV nim 's are positive reflecting the presence of 
symmetrical trions lying on the same lattice site, in the 
A^n, to) phases we find N n ^ m > but N miP < as well 
as N p _ n < reflecting the fact that the atoms of specy p 
lie on neigboring sites where the pairs [n, to) sit. In addi- 
tion, there remains to discuss the effect of the three-body 
losses |ll| which will reduce the lifetime of the trionic Ao 
phase, but are expected to have little effect on the or 
BCS phases. Therefore, provided that the temperature is 
low enough, current available experiments could achieve 
a BCS pairing instability in the (2,3) channel at small 
density or a A 7T (2, 3) phase for larger densities. 
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